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Equations of Motion of Systems with 
Internal Angular Momentum 

Manuel Dorado* 

CENITA, S.A., Miguel Yuste, 12, 28037 Madrid, Spain 
Abstract 

Using the Euler's equations and the Hamiltonian formulation, an attempt has been made to obtain the equations 
of motion of systems with internal angular momentum that are moving with respect to a reference frame 
when subjected to an interaction. This interaction involves the application of a torque that is permanently 
perpendicular to the internal angular momentum vector. 

mass) L, whose centre of mass moves at a constant 
velocity v with respect to a reference frame which 
can be defined as soon as the interaction in the 
cylinder starts. It is aimed to obtain the equations 
of motion to describe the dynamical behaviour of 
the system from the instant that it undergoes an 
interaction, by applying a torque M that is perma- 
nently perpendicular to the internal angular mo- 
mentum (See Fig.l). 

When solving the problem, the following points 
are taken into account: 



INTRODUCTION 

The angular momentum of a many-particle sys- 
tem with respect to its centre of mass is known as 
the "internal angular momentum" and is a prop- 
erty of the system that is independent of the ob- 
server. Internal angular momentum is therefore 
and attribute that characterizes a system in the 
same way as its mass or charge. In the case of a 
rigid body and particularly in the case of an ele- 
mentary particle, the internal angular momentum 
is also referred to as " spin" . 

To the author's knowledge, the dynamical be- 
haviour of systems with internal angular momen- 
tum has not been systematically studied within the 
framework of classical mechanics (See References 
[1]). This paper approaches this type of problem 
through a model comprising a rotating cylinder, 
with constant velocity of rotation, around its longi- 
tudinal axis. The equations of motion are obtained 
using the Euler's equations and by the Hamiltonian 
procedure. 

Results coincide when the problem is solved us- 
ing vectorial algebra or Lagrangian formalism (un- 
published). However, as a result of changes in "per- 
spective" , each method uncovers new peculiarities 
regarding the intimate nature of the system's be- 
haviour. 

FORMULATION OF THE PROBLEM 

Let us consider a rigid cylindrical solid with internal 
angular momentum (with respect to its centre of 



(a) The cylinder will continue spinning about its 
longitudinal axis at a constant angular veloc- 
ity uj throughout all the movement. In other 
words, its internal angular momentum mod- 
ule is constant. The energy that the cylinder 
possesses as a result of its rotation about its 
longitudinal axis is consequently considered to 
be internal and does not interfere with its dy- 
namical behaviour. 

(b) The derivative of the internal angular momen- 
tum, L, with respect to a frame of axes of iner- 
tial reference (X, Y, Z) satisfies the equation 




XYZ 




+ nxL (i) 



X'Y'Z' 
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in which f2 is the rotation velocity of the frame 
linked to the solid (X' , Y' , Z 1 ) about the frame 
of inertial reference axes (X, Y, Z). 

(c) Only infinitesimal motions are considered that 
are compatible with the system's configura- 
tional variations brought about by the applied 
torque. 



(d) The total energy of the system is not explicitly 
dependent on time. 

(e) No term is included that refers to the potential 
energy. According to hypothesis, the applied 
torque is a null force acting on the system and 
the possibility of including a potential from 
which the applied torque derives is unknown. 




Figure 1. Formulation of the problem 

HOW TO APPROACH THESE PROB- 
LEMS THROUGH THE EULER'S EQUA- 
TIONS? 

The figure 1 shows a cylindrical rigid body with 
angular momentum, L, about its longest axis. We 
should consider two referential frames: one associ- 
ated with the solid (X', Y' , Z 1 ), having the Z 1 axis 
on the direction of the body's longest axis. We will 
refer this system of coordinates to a inertial frame 
(X, Y, Z). 

In a given instant, a torque M is applied on the 
rigid body. 

We call f2 to the rotational velocity of the sys- 
tem of coordinates associated with the solid (X 1 , 
Y', Z'), viewed from the inertial frame (X, Y, Z). 
From here, we will consider the system of coordi- 
nates (X 1 , Y', Z 1 ) all along the work. 

Assuming this situation, it is known that 




XYZ 




+ flxL (2) 



X'Y'Z' 



It is clear that we have two different velocities 
acting: 



(a) uj: the rotational velocity of the body around 
its longest axis (let us call it the intrinsic ro- 
tational velocity of the body). 

(b) Q,: the rotational velocity of the system of co- 
ordinates associated with the body (X' , Y' , 
Z') viewed from the frame of inertia (X, Y, 
Z). 

We can write L and M, refcred to (X 1 , Y', Z 1 ), 
as follows: 



L = huiei + I 2 tu 2 e 2 + / 3 w 3 e*3 (3) 
Q = Qiei + VL 2 e 2 + Q 3 e 3 (4) 
M = Miei + M 2 e 2 + M 3 e 3 (5) 

So, the equation (2) takes the expression 

+ I 3 ^3^2 - hu 2 &3 = Mi 
I 2 Cj 2 + hujitlz - I 3 w 3 ^i = M 2 (6) 
I 3 u} 3 + I 2 lu 2 Q,i - I1LO1Q2 = M 3 

These equations are known as modified Euler's 
equations and can also be found in (l.e). 

In the case treated above, we must substitute the 
angular momentum 



in (6): 



L = I 3 w 3 e 3 



73^3^2 = Mi 
/ 3 W3^i = M 2 

I 3 uj 3 = M 3 



(7) 



(8) 



And then we get 



2 = 
fti = 
lo 3 = 



Mi 

I3U3 
M 2 

I3U3 
M 3 



(9) 



If the torque is on the axis Y', Mi — 0, M 3 = 
and: 



^2 = 



ill 

lo 3 



M 2 




(10) 



Where Sli coincides with the known velocity of 
precession. 

We have obtained that the body moves along a 
circular trajectory with a rotational velocity fii, 
but this is not enough to determinate the radius of 
the orbit. 

The principle of conservation for the energy 
should lead us to the expression for this radius. 

Assuming our initial hypothesis, the forces are 
applied perpendicular to the plane of movement of 
the solid and so, these forces do not produce any 
work while the cylinder is moving. 

We have two conditions that have to be satisfied 
by the body while rotating: 

(a) n = M. 

(b) The total energy should remain constant. 

Before the torque is applied the energy of the 
cylinder takes the following expression: 



E : 



1 2 1 r 2 

-mv Q + -Iuj 



(11) 



When the torque is acting on the cylinder about 
the line defined by the unit vector e 2 , it rotates 
along a trajectory, generally defined by r(t) and 
9{t) that, in this case, coincides with and its 
energy takes the following expression: 



(12) 



In this particular case our axes are principal axes 
of inertia, therefore we have ^gf- — Mi , and then, 
j- t {mr 2 6) = 0. 

And therefore mr 2 8 = const. 

In the particular case in which 9 = ft is a con- 
stant, r has to be a constant (that is, f has to be 
zero) in order to satisfy this equation. We can con- 
clude that the body moves in a circular orbit. The 
energy takes the following expression: 



E = -mr 2 n 2 + -Iu 2 



(13) 



Both expressions (11) and (13) have to represent 
the same energy. Therefore, comparing the two ex- 
pressions we can obtain the radius of the orbit: 



n 



(14) 



Looking at the results arising from the previous 
study, we can say that the body rotates following 
a circular orbit, with a radius given by the last 
expression. This circular trajectory implies that 
the velocity vector must precess jointly with the 
intrinsic angular momentum vector. From this last 
conclusion it can be proved that the force needed 
to cause the particle to draw a circular trajectory 
is expressed as: 



F = mP x f2 



(15) 



When the internal angular moment of the system 
does not coincide with one of the principal axes of 
inertia, the treatment of the problem is much more 
complex, but as it will be seen in the next section, 
the results given above area completely general. 

HAMILTONIAN FORMULATION. 
EQUATIONS OF MOTION AND THEIR 
SOLUTION 

The equations of motion are obtained by the Hamil- 
tonian formulation, where the independent vari- 
ables are the generalized coordinates and moments. 
To do this, a basis change of the frame (q, q, t) to 
(p, q, t) is made using the Legendre transformation. 
From the function 



a system of 2n + 1 equations is obtained. 



(16) 



dH 

_dL 

dt 



-Pi = 

dH 
~dt 



dqi ' 



(17) 



If the last equation in (17) is excluded, a system 
of 2n first order equations is obtained, known as 
canonical Hamilton equations. 

By substituting H, the first order equations of 
motion, are obtained. 

It should be noted that the Hamiltonian formula- 
tion is developed for holonomous systems and forces 
derived from a potential that depends on the po- 
sition or from generalized potentials. A torque is 
applied to the present system. The result of the 
forces is zero on the centre of mass and, therefore, 
it is meaningless to refer to potential energy. 



On the other hand, the Hamiltonian function The term — ^ normally appears in this radial 
concept does remain meaningful. equation of motion and represents the force derived 

By using polar coordinate in the movement from a potential. In other words, 
plane, 

* = v °= r§ ( 18 ) mf = mr^-? (27) 

or 



The Lagrangian is 



1 . In the present case, this term does not exist. 



L ~ 2 m(r +r ^ ( 19 ) To complete one's knowledge of the system's de- 

zed moments are 
P r = mr; Pg = mr 2 6 (20) 



vclopment, those relations derived from the con- 
Thc generalized moments are , . , , , , , 

° stramts, have to be resorted to. 



so that 



m mr* 



1. E = imv 2 = const. => \v\ = const 

^ ^ (ti«) X yz ( d * ) X'Y'Z' 



The Hamiltonian H is introduced using the equa- 
t j on where the applied external torque 

H{p,q,t) = ^2q i p i -L(q,q,t) (22) ^ 

M= f ) (28) 
resulting m \ al 



XYZ 



H = P r r + P e 9 - (-mr 2 + -mr 2 6 2 \ by hypothesis it is known that 

^ 2 2 ' / / = and from (1) it is concluded that 

If the generalized velocities are substituted by ( dZ\ - O x L 

the generalized moments, the following is obtained V dt ) x'Y'Z' 

and it is obtained that O = 

p2 p2 L 

jj — r _|_ 8 z'23) This vector represents the rotation velocity of the 

2m 2mr2 frame of axes linked to the cylinder (X', Y', Z') 

The first pair of Hamilton equations is about the frame of inertial references axes (X, Y, 

Z). 

f = — = Q = — = ^ In polar coordinates, the variable defining the 

®Pr m dPe mr rotation about the frame of axes is 9 and it is con- 

The second pair of equations is eluded that 

M 

_ Pr = dH = _p$_ _ Pe = ^L = {) e = n = T (29) 

dr mr^ 00 

substituting in the radial equation of motion, re- 
The second of these equations shows that the an- suits that 
gular momentum J is conserved. m ~ _ mr Qi _ mr Q2 ^q-j 

Pg = J = const. (24) 

Moreover, Pg=const. and at the initial moment 
The first gives the radical equation of motion, equals 

Pg — mr 2 6 = mrv (31) 



J2 

P r =mr = - (25) 



as v is constant, it is concluded that r is also 



as J = Pg = mr 2 6\ by substituting it results that constan t- Fin ally it is found that 



P r = r ^=mre 2 (26) r = S and ™ = ™* ^ 



Further considerations 

The Hamiltonian is independent of 6. This is an 
expression of the system's rotation symmetry or, 
in other words, there is no preferred alignment in 
the plane. 

The equation Tjf- = means that the energy of 
the system remains invariable if it is turned to a 
new position without changing r, P r or Pg. This 
aspect is clearly proved. 

Worthy of mention in this case is the fact that 
the Hamiltonian equations do not provide, in this 
type of problem, complete information on the ra- 
dial movement. The reason for this is that the cen- 
tral force is a function of 6 which is a coordinate 
that does not appear in the Hamiltonian (it can 
be ignored) and can only be calculated using the 
constraints. 

This proposal is valid for any property of the 
particle that displays vectorial characteristics and 
is sensitive to an external torque of the type de- 
scribed, without the need to ascertain the real phys- 
ical nature of that property. The proposal can be 
extended to equivalent problems, despite their not 
involving the presence of internal angular momen- 
tum as a central characteristic. 

Poisson brackets. The Poisson bracket [Pg,H] is 



[Pe,H} = 



dP e dH dP e dH 



86 8Pg dPg 86 



dH 
~86 



as Qjg- = 0, it is concluded that [Pg, H] = 0. 

CHARGED PARTICLE IN A MAGNETIC 
FIELD 

Let us apply this theory to the specific case of a 
charged particle with spin s and magnetic momen- 
tum fl, moving at a velocity i? in a uniform mag- 
netic field of flux density B. 

As it's well known, the interaction between B 
and fl makes a torque to act upon the particle. This 
situation verifies every condition we have defined in 
our hypothesis, and so, the behaviour of the particle 
must be in agreement with the theoretical model we 
have developed above. 

Considering the most general of the cases, let us 
suppose that fl forms unknown angles 6 with B 
and a with vq. This situation is represented in the 
following figure (See Fig. 2). 




Figure 2. Charged particle in a magnetic field. 

It is clear that the rotation of Uq causes rotation 
of the i?o component on the XY plane, v§ xy . Fur- 
thermore, according to the theory here presented, 
the particle will trace a helix and the Vq compo- 
nent in the direction of the Z axis, v$ z , will not be 
affected. 

In this case the angular momentum L = s. 

If the particle stops inside B, the rotational 
velocity of the angular momentum would be ex- 
pressed by: 



n 



\jixB\ 



\M\ 
\L 



^ (33) 



Obviously, if the angle a formed by s and v 
remains constant, the angle 7 formed by s xy and 
voxy will also remain constant. The velocity of 
can be calculated merely by calculating the velocity 

= ssin6> (34) 



?xy 



and 



Q>V0xy — Qsqx 



\fl\\B\smd 
L sin 6 



I^B 



(35) 



We can conclude from this expression that the 
angular velocity of the rotation of the particle de- 
pends on the magnetic momentum jl, the magnetic 
field B in which the particle is immersed and the 
spin s of the particle, and it is independent of the 
relative positions of fl, B and v$ . 

In the case in question, the angular velocity of 
the rotation of the particle in B is expressed by: 



ft = 7 B 

'2m 



(36) 




According to the theory developed above, the 
particle will draw a circular trajectory. We can 
determine the radius of the orbit drawn by the elec- 
tron under the influence of a magnetic field 



jeB 



(37) 



In the specific case in which the particle is an 
electron, 7 has a value of 2 and the radius of the 
orbit is calculated from the following equation: 



^0 
O 



~eB 



(38) 



From the conclusions summarized above, we 
state that the particle is submitted to a central 
force affecting the charged particle inside the mag- 
netic field. To calculate its value, we only need to 
recall the formula: 



F = mpQ x Q 



(39) 



and substitute f2 for its value previously obtained. 
Then we get: 



F = mv§ x 7 B 

2 m 



(40) 



Once again, if we apply this to the case in which 
the particle is an electron, the gyromagnetic factor 
is 2, and hence 



F — cvq x B 



(41) 



The behaviour predicted by the theory here 
shown of a charged spinning particle penetrating 
into a magnetic field with velocity ^, matches the 
one that can be observed in the laboratory, and 
the force to which it will be submitted is the well 
known Lorentz force. 

EXPERIMENTS CARRIED OUT 

Airmodel with spinning disc 

The airmodel is provided of a strong angular mo- 
mentum by mean of a high speed spinning disc at- 
tached to it by exerting different torques on it, with 
the built in flight systems, we achieve modifications 
in its flight trajectory. 

The remote control flying model was designed 
specifically to verify this theory (See Fig. 3). 




Figure 3. Airmodel with spinning disc 

Basically it consists of a disk that can rotate, 
attached on top of the flying model, which provides 
the necessary support and speed. 

Manoeuvring the flying system, aileron to turn 
right and left, elevators to shift upwards and down- 
wards and rudder, we can obtain the adequate 
torques to affect this motion. 

The flying model is controlled by means of a radio 
equipment affecting the normal moving elements 
and the rotation of the disk. 

The fist flight was performed without the disk to 
confirm the normal evolution of the model. 

Once the normal behaviour of the model in this 
conditions was confirmed, the disk was attached on 
top of the model. 

With the disk still the system behaves in a sim- 
ilar way that it did without the disk, although the 
stability was slightly affected. 

With the disk rotation pointer clock way when we 
actuate the ailerons to the left, instead of drawing 
a horizontal circumference, it actually goes up even 
to the point of performing a loop. 

Actuating the elevators upwards causes the 
model to turn right, when we actuate the elevators 
downwards the model turns to the left. 



Spinning top with double suspension 
assembly 

In this experiment it is proved that the spinning 
top cam precess with a non zero radius. 

It can be observed as well its extraordinary sen- 
sitivity to the torques applied on it, when these 
torques have the same direction of the vertical axis, 
by modifying the angle between the symmetry axis 
of the spinning disc and the earths plane. 

For this experiment, a top is mounted with dou- 
ble suspension (See Fig. 4) whereby it can behave 
according to the general theory or adopt the type 
of motion predicted by the classical treatment. 




Figure 4. Top mounted with double suspension. 

The top is allowed four possible degrees of rota- 
tional freedom: 

a) AA' axis: The top rotates about its axis of 
symmetry. 

b) BB' axis: The top will rotate about this axis 
if nutation is caused. 

c) CC axis: The top will rotate about this axis 
if there is precession. 

d) DD' axis: The system will rotate about this 
axis if it behaves according to this theory when 
we fixed the CC axis to the support. 

The gravitational torque follows the direction of 
the BB' axis and it produces an increase in angular 
momentum AL which is vectorially added to the 
angular momentum L of the top, which is in the 
direction of the AA' axis. 



Logically, the top and both the AA' and BB' 
axes could occupy any position with respect to the 
support frame. 

We provide the top with a constant slow velocity, 
Vq. Due the top being mounted on the support 
frame, the system will rotate about the DD' axis. 
As Vq is slow, the centripetal acceleration due to 
rotation can be considered negligible. 

According to the classical explanation, the top 
should precess and even achieve a nutation motion 
on the support, with an independent slow rotation 
motion of the system about the DD' axis. 

The general theory establishes that the velocity 
of the top, i>o, the orbital radius and the preces- 
sion velocity of the top are related by the following 
expression: 

Surprisingly this means that if f2 and R are con- 
stant (the support frame is rigid) there would be 
only one velocity v$ for the movement of the top 
and a single rotation velocity for the top-support 
system about the DD' axis, independently of the 
impulse given to the support in an attempt to 
achieve the desired velocity. 

For our experiment, we create a rotation veloc- 
ity for the top about its axis of symmetry, this axis 
(AA') being maintained initially in a horizontal po- 
sition although it can really be in any position. 

When the top has reached a high rotation veloc- 
ity about its axis of symmetry, the support is given 
an impulse so that the top reaches a velocity V and 
the system is then left to evolve freely. For now 
on, it is also under the influence of a gravitational 
torque. 

If the system behaves according to classical the- 
ory, precession of the top should not occur with 
respect to the support. 

If on the other hand its behaviour is that de- 
scribed by this theory, precession should occur with 
respect to the DD' axis. 

The results of the experiment confirm that the 
top, in agreement with the theory, move about the 
DD' axis with a velocity 

= R x ^0 

and is independent of the impulse it may have re- 
ceived. 

In fact the following situations can occur: 



a) That the impulse received may be bigger than 
that required to reach a velocity v . The ve- 
locity of movement of the top will in all cases 
be i7q. However, the angle formed by the axis 
of symmetry (A A') and the vertical axis (CC) 
decreases and there is an increase in potential 
energy 

It should be remembered that f2 is constant and 
independent of the angle formed by the AA' and 
CC axes while L remains constant. 

b) That the impulse received is equal to that re- 
quired to reach a velocity u . In this case the 
position of the top does not vary with respect 
to the support. The system rotates about the 
DD' axis. 

c) That the impulse received is less than that re- 
quired to reach a velocity Vq. In this case the 
system continues to move with a velocity Uq, 
but the angle formed by the AA 1 and the CC 
axes increases. 

Pendulum with internal angular mo- 
mentum 

Two experiments are performed with this device: 

a) Measure the curvature of the trajectory at the 
very first instants of its evolution, where we 
assess that it is close to a circular trajectory. 

b) Check that the final trajectory is an ellipti- 
cal one in which the perihelio of the orbit ad- 
vances. 

This experiment is performed with a pendulum 
that consist of a fibre-glass sphere. The sphere 
houses a disk that can rotate around its transversal 
axis (See Fig. 5). 

The pendulum is suspended from a wire attached 
to a point one centimeter distant of the edge of the 
disk shaft. 

If we observe the evolution of the pendulum when 
the inner disk stands still, it is the ordinary one. 

We force the disk to rotate until the maximum 
speed is reach. At that moment we let the sphere 
free. The system is under the influence of a torque 
caused by its weight and the wire strain. If we 
measure the trajectory at the very first instants of 
its evolution, it happens to be a round one. 



Letting the sphere to move, after while it follows 
an elliptical orbit, on which the pcrihclium shifts in 
time, in good agreement with this theory (See Fig. 
6). 

This phenomenon is similar to the one occurring 
in the shift of a planet orbit perihclium. And it 
is similar, as well, to the Larmor precession of a 
charge in the presence of a magnetic field. 




Figure 5. Pendulum with internal angular momentum. 

As in the previous experiment, if we reverse the 
direction of the disk rotation the orbit pcrihclium 
will also shift in the opposite way that it did before. 




Figure 6. Pcrihclium shift 

Theoretical explanation 

— = M (42) 
dt v ; 

M = fxF (43) 



In agreement with this theory, the pendulum is 
affected by the force 



F = mv x 0, 



The torque will be 



M = rx mv x Q, = L x Cl 



(44) 



(45) 



(It can be proved that the associative property of 
the vectorial product is satisfied in this case) 
Substituting (45) in (42) we obtain: 



dL r n 

dt 



(46) 



This one is the equation of the motion of a vector 
with constant module \L\, that precesses around 
the axis defined by ft, with angular velocity |f2|. 

Larmor Effect (Samples under the influence of mag- 
netic fields) 

In agreement with this theory an electron em- 
bedded inside an atom in the presence of a mag- 
netic field B, is submitted to the central force 



F = mv x f2 
which adopts the expression 



(47) 



F = mvx(i = mvx — = -vxB (48) 
mc c 

and is known as the Lorentz Force. 

As the most general expression, we can write: 



dL 
dt 



f x mv x Q = L x 



(49) 



(It can be proved that the associative property of 
the vectorial product is satisfied in this case) 

This is the equation for the motion of L, which 
rotates around the vector B with an angular veloc- 
ity: 

q 



-B 



mc 



(50) 



The general expression of the central force, which 
includes the Lorentz force as a particular case, al- 
lows us to explain the Larmor effect and its mean- 
ing in an easy and accurate way within the range 
of a wider and more important phenomenon. 



Magnetic spinning top 

With this experiment it is proved that the direc- 
tion way of the precession can be manipulated in 
accordance with the torque applied. 

For the this experiment we have a magnetized 
rod, a aluminium cone and a magnet. 

We attach the rod to the cone. We have build 
up a magnetized spinning top (See Fig. 7). 



B 




Figure 7. Magnetic spinning top. 

We observe the behaviour of the magnetized 
spinning top in the presence of a magnetic field. 
The interaction between both fields produces a 
torque which affects the spinning top making him 
to draw a circular trajectory. 

The spinning top will not collapse on the mag- 
neto as long as it has an angular momentum. 

EXPERIMENT SUGGESTED TO TEST 
THIS THEORY: INTERACTION OF A 
HOMOGENEOUS MAGNETIC FIELD 
AND A PERPENDICULAR GYRATING 
MAGNETIC FIELD, WITH A PARTICLE 
WITH SPIN AND MAGNETIC MOMENT 

Formulation of the problem 

The particle with spin, S, moves with velocity v. It 
first passes into a homogeneous magnetic field B' , 
inducing the magnetic moment, jx, of the particle, 
polarized in the direction of the field. It later enters 
a region in which a homogeneous field, B, and a 
gyrating field, Hq, with rotation frequency, 0, are 
superposed. Fields B' and B are parallel and in 
fact could even be the same field, although their 
function is different in each region. 



Reference frame (See Fig. 8). 
XYZ Fixed system in the particle but which 
does not rotate with this. 

B is on the Z axis. 
Ho is on the XY plane. 
X'Y'Z' Fixed frame in the particle which pre- 
cesses with this. 

S and p, are on the Z' axis. 
V is on the Y' axis. 
XqYqZq Reference frame from which observa- 
tions are made. 

B is on the Zq axis. 

Hq is on the X Y plane. 




Y 



Figure 8. Reference frames. 

The rotation velocity, 0, of the frame of axes 
linked to the solid, X'Y'Z', will be expressed within 
this system as: 

Q = fii? + Q 2 j' + Q 3 k' (51) 

This rotation velocity, O, can also be expressed 
with respect to the frame of reference axes, XYZ 
by using Euler angles. 

Qi = <j> sin sin tp + cos \p 

Q.2 = sin 6? sin — ^ cos V' (52) 

3 = <fi cos + ip 

Our particle can also rotate with respect to the 
frame of axes, X'Y'Z', and the most general ex- 
pression to describe this angular velocity is: 

Q = uix'i + WY'j' + ojz'k' (53) 



It is proved that ip = u>z' ■ The proof is easy, for if 
the Z' axis is fixed, i.e. both <p> and are constant, 
the rotations about the Z' axis are rotations of the 
X' and Y' axes in the X'Y' plane. Hence, the Euler 
angle that reports the rotation is ip. It can then be 
concluded that 

ip = u) Z ' (54) 

By agreement, but without loss of generality, let 
us consider that the particle has an angular mo- 
mentum, S, that is constant in module. 

Frame of axes linked to the solid 

In the definition of the frame of axes linked to the 
solid, we can either choose: 

a) The frame of axes strictly accompanies the 
particle in the rotation which, in classical 
terms, gives its internal angular momentum. 

b) The frame of axes linked to the solid is de- 
fined by a parallel axis to the internal angu- 
lar momentum of the particle, and the other 
two are perpendicular to each other and are 
in a plane which is perpendicular to the angu- 
lar momentum. Internal angular momentum 
is understood to be a quality of the particle 
whose mathematical characteristics coincide 
with those of the angular momentum, without 
considering the physical nature of this quality. 
In other words, no hypothesis is formed as to 
whether or not the angular momentum implies 
rotations of the particle under study. 

For the purposes of our problem, definition (b) 
has been used. 

However, it should be remembered that our par- 
ticle can reach a rotation velocity, u, as a result of 
the interactions to which it can be subjected. 

Analysis of the interactions 

Given the nature of our problem, the particle pen- 
etrates a homogeneous magnetic field, B, inducing 
the magnetic moment, jl, of the particle, polarized 
in the direction of the field. 

Hence, when the particle penetrates the magnetic 
field, B (parallel to B' , and can even coincide with 
it) and Hq, the magnetic moment, /2, is only (ini- 
tially) sensitive to the gyrating magnetic field, H . 



This, consequently, will be the first interaction that 
we analyze. 



Interaction with the magnetic field, 

H 

The interaction of the gyrating magnetic field, H , 
with the magnetic moment, ft, does not modify the 
energy of the system, as both are permanently per- 
pendicular. 

The energy of the interaction is expressed 

E H = H - n = Hon cos - = (55) 




T 



Figure 9. Interaction with the magnetic field, Hq. 



The interaction between Ho and [2, is equivalent 
to the action of a torque, Th, on the particle, where 



T H = fix H 



(56) 



0, 



S 






(59) 



Th is perpendicular to the X' and Z' axes, and 
is located on the Y' axis. 

By using the previously obtained Euler equations 
(9) and applying them to our problem: 



UJ Z ' 
^2 



Ms 
h 

s 

Mi 



= = -0 
S 

= 



(57) 



so that ip = in all cases, given that M 3 = at all 
times. We can avoid tp by making it equal to zero 
in all case (See Fig. 9). 

The Euler equations depending on the Euler an- 
gles (52), are simplified with the result that 



n 2 = 
n 3 = 



5* 

) sin 6> = 

)COS0 



(58) 



in other words, the rotation of the frame linked to 
the solid on X'Y'Z', Oi, is constant and equal to 6, 
where 9 is the rotation of the frame of axes linked 
to the solid with regard to XYZ expressed in the 
Euler angles. 

We have obtained the rotation velocity, 9, but we 
have to calculate the rotation radio r(t) if we want 
to know the trajectory. To do this we will use the 
following equality between differential operators, 



d_ 
~dt 



dt 



where 

4t is the derivative in the inertial frame, in our case 
XqYqZo, 

^ is the derivative in the non-inertial frame, in our 
case X'Y'Z', 

(l is the angular rotation velocity of the frame of 
axes linked to the solid in relation to the inertial 
frame of axes. 

We will apply this equality to r, a vector of posi- 
tion of one frame with respect to another, and for 
convenience we will express all the vectors in the 
X'Y'Z' frame. 



from which it is deduced that 



d*f dr -> 

— — = fi x r 

dt dt 



(60) 



If we represent the equation in components, we Interaction with the magnetic field B 

have 

At the moment the circular trajectory is com- 
menced, 9 is no longer null and the interaction with 
the field, B, begins (See Fig. 11). 



rx' + r Z '^2 - r Y '^3 = v x > 
rz 1 + ryfii - rx'^2 = v z < 



(61) 



Taking into account the initial conditions (See 
Fig. 10), n 2 = and fi 3 = 0. 



X X' 
Line of nodes 



r = 
n 2 = 



fl x B 
fj,B sin 9 



(63) 
(64) 




Figure 10. Initial conditions. 

The interaction does not modify the energy of the 
system, as explained beforehand, so the module of 
v has to remain constant. 

\V\ =v = constant. 

and, vx< = 0, vyi = i>, vz* = 0. 

As the energy of the system cannot vary, fx 1 = 
fy = fz 1 =0. If they differed from zero tangen- 
tial accelerations would be involved and therefore 
a variation in the total energy of the system. 

These conditions simplify the equations of our 
system (61) and we obtain 



rz> = - 



(62) 



This is to say, the particle and the frame of axes 
linked to the particle describe a circular trajectory 
with respect to the observation frame. The move- 
ment plane is perpendicular to H n . 



Figure 11. Interaction with the magnetic field. 

M 3 continues to be null and, therefor, so do tp 
and tp, which do not vary, but Q 2 is now different 
from zero. 



fii = 6 

Q.2 — <j)sm8 = 

Cls = (j) cos 9 



fiB sin 8 



(65) 



We finally obtain, 



O, = 



S 
fj,B 

~T 

[iB cos 9 



(66) 



These are the expressions of the rotation of the 
frame of axes linked to the solid in the base that 
defines the trihedral linked to the solid. 

We can also express this same rotation velocity of 
the frame of axes linked to the solid with respect to 
the reference frame, by means of the Euler angles, 



but in the base defining the trihedral of reference 
frame. 



ujx = cos <p + tp sin 9 sin < 
uj y — sin <f> — ip sin 9 cos < 
loz = ip cos 9 + (j) 

As V> = V = 0, we obtain 



(67) 



UJ X 
UJy 

LO Z 



= V COS ( 

= 9 sin c 



(68) 



These expressions clearly reflect the behaviour of 
the frame of axes linked to the solid and of the solid, 
with respect to the reference axes. 

Visualization of the movement is assisted by con- 
sidering the following specific situations = 



i) For 







uj x 

UJy 
UJ Z 



We would have a rotation about the X axis, 
j) For 4> = f 

ujx = 



UJy 
UJ Z 







We would have a rotation about the Y axis. 
The rotation velocity in the XY plane is 

wxy = 9(cos 9i + sin 9 j) 



(69) 



\uxy\ — 



Therefore, the total movement is a circular tra- 
jectory of radius, rz>, perpendicular to the field, 
Ho, which in turn spins about the field, B, with 
velocity 4> (classically know as Larmor rotation fre- 
quency). This explains why the frequency of the 
oscillating field, Ho, has to coincide with the Lar- 
mor frequency. 



The behaviour of the particle indicates that it is 
subjected to two central forces 



F\ = mv x Oi = mv x 9i' 



(70) 



mv x 4>k (71) 

From an energy point of view, we can confirm 
that the magnetic field, H , does not modify the 
energy of the system, and that the magnetic field, 
B, cyclically modifies the potential energy of the 
particle, since 



Eb = B ■ flcost 



(72) 



with 



0(t) = 



but with its kinetic energy unvaried. 

Additional considerations 

The reader will have already noted the similarity 
between this behaviour and the Larmor effect. 

It is know that this effect is produced by applying 
a magnetic field, B, to a particle of charge, q, mov- 
ing in an orbit around a fixed specific charge, q' . 
The result is a precession of the trajectory around 
the direction of the applied magnetic field, with a 
precession velocity, ujl = know as Larmor fre- 
quency, with the proviso that the cyclotronic fre- 
quency is directly obtained by this procedure. 

This statement can be expressed by saying that 
the angular momentum vector of the particle with 
respect to the rotation axis, referred to as orbital 
angular momentum, precesses about the direction 
of the magnetic field, B. 

From our procedure, it is particularly easy to 
reach the same conclusions by considering that 



at 



L = fx mv 



(73) 



By substituting our development, F = mv x fi, 
in (73) we obtain 



— = r x (mv x ft) 
dt v ' 



(r x mv) x (} = L x f2 



(74) 



Z 



B Homogeneous 
Magnetic Field 
and 



B" InhoniDgeneous Magnetic Field 




B' Inhomogeneous Magnetic Field 



Hq Oscillating 
Magnetic Field 



Figure 12. Rabi's experiment. 
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Figure 13. Experiment suggested. 



(the triple vector product does not generally 
fulfill the associative property, but it can be de- 
mostrated in our.) 

Equation (74) describes the angular momentum 
precession about the magnetic field with a preces- 
sion rate f2. 

The Larmor effect results from subjecting the 
particle to the Lorentz central force, because it is 
within B. From the following, the Lorentz force is 
equivalent to our force, given that 

F = mvCl = mv^- = mv—B = evB (75) 

Which demostrates that both movements are 
equivalent. 



Experimental conditions 

This experiment can be carried out by emulating 
the rotating magnetic field using a radio frequency 
magnetic field, just as Rabi did in his experiments. 

Rabi's experiment (see References (2)) consists 
of a collimated particle beam that crosses an inho- 
mogeneous magnetic field. It later passes through a 
region where a homogeneous and a radio frequency 
magnetic field are superposed, and finally passes 
through an inhomogeneous field that refocuses the 
beam towards the detector. 

The inhomogeneous fields separate the beam into 
different beams according to their magnetic mo- 
ment (the dependence this on the spin) as in an 



experiment of Stern-Gerlach. When leaving the 
first field, these beams later pass into the second 
inhomogeneous field, which refocuses the beams to 
the detector. By adjusting the second inhomoge- 
neous magnetic field we will obtain refocusing con- 
ditions for a beam or group of beams. 

In the central part of the experimental arrange- 
ment, the homogeneous magnetic field is super- 
posed on the radio frequency field. In this region, 
the spin and magnetic moment are deflected with 
respect to the constant field when the oscillating 
field frequency approaches the Larmor precession 
frequency. This process is know as nuclear mag- 
netic resonance. 

After deflection, the atom is on another level 
which will not fulfill the refocusing condition and 
a decrease in the beam intensity will be observed 
in the detector. This procedure is used to study 
nuclear spin, nuclear magnetic moments and hy- 
perfine structures. However, the behaviour of the 
particles in Rabi's experiment is theoretically jus- 
tified in our exposition and it is therefore unnec- 
essary to use inhomogeneous magnetic fields (as is 
required in Rabi's experiments). 

Our experiment only requires the presence of a 
homogeneous magnetic field to induce the magnetic 
moment, jl, of the particle, polarized in the direc- 
tion of the field. 

In such a way that: 

1° Only those particles whose Larmor frequency, 
<t>, coincides with the radio frequency will 
abandon the rectilinear trajectory. 

2° The trajectory described by these particles 
should coincide with that of the theoretical so- 
lution given in this article while they remain 
within the region in which the two fields co- 
exist. According to this solution, the parti- 
cle abandons the rectilinear path and the XY 
plane. 



CONCLUSIONS 

When a system with internal angular momentum, 
moving at constant velocity with respect to a ref- 
erence frame, is subjected to an interaction of the 
type under study, i.e., a torque perpendicular to the 
internal angular momentum vector, it will begin to 
trace a circular path of radius 

v 

r= ri 

where = M. 

Its behaviour is equivalent to that produced by 
subjecting the cylinder to a central force 

F = mv x (l 

As it can be seen in figure 14, the trajectory 
which particle follows is the trajectory II, while the 

trajectory I is the intuitive one, but as it has been 
shown in this paper it is not the real behaviour of 
the particle. 
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3° When the particle abandons the region of co- 
existing fields, it will follow a rectilinear move- 
ment but will not reach the detector. 

4° This procedure facilitates selection of the dif- 
ferent angular moments because of their de- 
pendence on the Larmor precession. 




Figure 14. Classical trajectory and trajectory predicted in this article. 
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